Atomic Motion in Single H 2 and D 2 Molecule Junction Induced by Phonon Excitation 
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We have investigated Au atomic contacts in Hb and D2 environment by conductance measurement 
and dl /dV spectroscopy. A single H2 or D2 molecule was found to bridge Au electrodes. In the case 
of the Au/Hh/Au junction, symmetric peaks were observed in dl/dV spectra, while they were not 
observed for the AU/D2/AU junction. The shape of the peaks in dl /dV spectra originated from the 
structural change of the single molecule junction induced by the phonon excitation. The structural 
change could occur only for the A11/H2/AU junction. The difference in the two single molecule 
junctions could be explained by larger zero point energy of AU-H2 vibration mode than that in the 
Au/Eb/Au junction. 

PACS numbers: PACS numbers: 73.63.Rt, 73.40.Cg, 73.40.Jn 
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INTRODUCTION 

Nano junctions using single molecules have attracted 
wide attention in recent year [HQ. In order to inves- 
tigate the atomic scale single molecular junctions, it is 
important to characterize a single molecule bridging be- 
tween metal electrodes using some spectroscopic mea- 
surements. As for the spectroscopy of the single molecule 
junctions, point contact spectroscopy (PCS) [3, lal , in- 
elastic tunneling spectroscopy (IETS) @, 0, H,M TU and 
action spectroscopy have been investigated for the single 
molecule junctions (llj . In these spectroscopies, electron 
current through the single molecule junction is measured 
as a function of bias voltage. In the IETS and PCS, 
spectrum is obtained as an inelastic scattering event of 
electrons with sufficient energy to excite certain vibration 
mode. The vibration energy is determined from the bias 
voltage, in which conductance of the single molecule junc- 
tion abruptly changes. The difference between PCS and 
IETS is operating conductance regime. IETS operates in 
low conductance regime (G <Cl Go, where Go —2e 2 /h), 
while PCS operates in conductance regime close to 1 Go 
i, H E3- The PCS and IETS have been applied for 
various single moleculeiunctions with H2, Cgo, benzene, 
alkane dithiol 0, i, H, (S S B S • 

Action spectroscopy of the single molecule junction 
was a recently developed technique to detect the vibra- 
tion modes of the single molecule junctions [ll[. The 
action spectroscopy observes the abrupt conductance 
change (peak in dl/dV spectra) originated from the 
structural change of the single molecule junction induced 
by the phonon excitation. The energy of the peak pro- 
vides the vibration energy of the single molecular junc- 
tion. The action spectroscopy of the single molecular 
junction was applied for the single molecular junctions 
with H2, O2, and CO ll|. In the case of single molecules 
adsorbed on metal surfaces, the action spectroscopy has 
some advantage compared to IETS [li , HI, 1|| . It can de- 



tect the vibration modes which could not be detected by 
IETS. In addition, action spectroscopy can be applied for 
the system, in which IETS signal is too weak to detect, 
since the signal of the action spectra is much larger than 
that of IETS. The action spectroscopy has been applied 
for various molecules on metal substrates. On the other 
hand, there is little study of the action spectroscopy of 
the single molecular junction, and thus, the character of 
this spectroscopy is not clear up to now. It is important 
to investigate the characters of the action spectroscopy 
of the single molecular junction. In the present study, 
the action spectroscopy has been performed for the Au 
atomic contact in H2 and D2 environment. This system 
has been already investigated by the conductance mea- 
surement, dl/dV spectra, and theoretical calculation by 



other groups [U El EI 



|20|] . However, IETS or PCS of 
this system has not been obtained up to now, and the en- 
ergy of the vibration mode was not fixed yet. The struc- 
ture of the Au atomic contact in H2 environment was not 
defined. In the present study, the atomic configuration of 
the Au atomic contact in H2 and D2 environment was in- 
vestigated by the conductance measurement, action spec- 
troscopy. Here it should be noticed that hydrogen has 
large zero point energy due to its small mass. Since the 
action spectroscopy detects the structural change of the 
single molecular junctions induced by the phonon exci- 
tation, it provides the information of the atomic motion 
of the single molecular junction. The effect of the zero 
point motion on the atomic scale structural change of the 
single molecule junction could be revealed by the action 
spectroscopy. In the present study, we have investigated 
the effect of the zero point motion on the atomic motion 
of the single molecular junction by investigation of the 
action spectroscopy of the Au atomic contacts in H2 and 
D2 environment. 
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EXPERIMENTAL 

The measurements have been performed using the me- 
chanically controllable break junction (MCB J) technique 
(see Ref. for a detailed description). A notched Au 
wire (0.1 mm in diameter, 10 mm in length) was fixed 
with epoxy adhesive (Stycast 2850FT) on top of a bend- 
ing beam and mounted in a three-point bending config- 
uration inside a vacuum chamber. In ultra high vacuum 
at 4 K, the Au wire was broken by mechanical bending of 
the substrate, and clean fracture surfaces were exposed. 
The bending could be relaxed to form atomic-sized con- 
tacts between the wire ends using a piezo element for fine 
adjustment. H2 or D2 gas was admitted to the contacts 
via a capillary. DC two-point voltage-biased conductance 
measurements were performed by applying a bias volt- 
age in the range from 10 to 300 mV. Every statistical 
data set was built from a large number (over 3000) of 
individual digitized conductance traces. AC voltage bias 
conductance measurements were performed using a stan- 
dard lock-in technique. The conductance was recorded 
for fixed contact configuration using an AC modulation 
of 1 mV amplitude and a frequency of 7.777 kHz, while 
slowly ramping the DC bias between -100 and +100 mV. 



RESULTS AND DISCUSSION 

Figure Q] shows typical conductance traces and his- 
tograms for Au contacts before and after introduction of 
H2. The stretch length was the displacement of the dis- 
tance between the stem parts of the Au electrodes which 
were fixed on the substrate. Before introduction of H2, 
conductance decreased in a stepwise fashion, with each 
step occurring at integer multiples of Go- The corre- 
sponding conductance histograms showed a peak near 
1 Go, which corresponded to a clean Au atomic con- 
tact. After introduction of H2, the conductance de- 
creased multi-stepwise fashion below 1 Go- The widths 
of the multiple steps were not well-defined, and showed 
various value less than 1 Go- The corresponding conduc- 
tance histogram showed a broad feature below 1 Go- The 
appearance of steps and features below 1 Go agreed with 
the previously reported results [l7j. Here, it should be 
noticed that the return conductance traces showed the 
steps around 0.1-0.3 Go, before making large contacts. 
The certain atomic configurations showing conductance 
values below 1 Go would be formed in both breaking and 
making contacts. 

The differential conductance (dl/dV) spectra were 
measured for the Au atomic contact in H2 environment 
at a conductance value of 0.1-0.5 Go- Three types of 
dl/dV spectra were observed, as shown in Fig. O The 
first "normal" spectrum (Fig. [21(a)) showed an increase 
in the differential conductance symmetrically around 45 
meV, and clear symmetric peaks were observed in the sec- 



ond derivative (dI 2 /dV 2 ). The second "peak" spectrum 
(Fig. [2Kb)) showed symmetrical peaks, and the third 
"non symmetric" spectrum (Fig.[2lc)) did not show any 
symmetric features. The third "non symmetric" spectra 
were most frequently observed for the Au contact in H2 
and D2 environment. The increase in conductance ob- 
served for the "normal" spectra could be explained by 
thephonon excitation at the single molecular junction 
B II 0, 0, 0, B 0, E El- The single-channel model 
predicts conductance enhancement below a transmission 
probability of 0.5 and suppression of conductance above 
this. In the present study, the spectra were measured 
for the junction having zero bias conductance below 0.5 
Go- The increase in the conductance agreed with the 
theoretical prediction. In the single molecular junction 
having low conductance, an additional tunneling chan- 
nel for electrons was opened, when the bias voltage was 
increased and crossed the threshold for excitation of a vi- 
bration mode. This opening of the new channel resulted 
in a sudden increase in the differential conductance at the 
threshold voltage. Here, it should be noticed that dl/dV 
spectra of the Au atomic contacts in H2 environment has 
not been obtained before 17, ID, ID, HI- I n the present 
study, we could measure the vibration spectra of this 
system by measuring the carefully prepared sample. The 
peak observed in the "peak" spectra could be explained 
by the abrupt switching between two slightly different lo- 
cal geometric configurations induced by the phonon ex- 
citation 11 1 . Action spectroscopy observes these peaks 



to determine the energy of the vibration mode. The non 
symmetric feature in the " non symmetric" spectra could 
be explained by interference of electron waves which scat- 
ter on defects or impurities close to the contact 0, 0]. 
Since the distribution of the defects or impurities and 
shape of the left and right electrodes would be different 
from each other, that is, non symmetric, the non sym- 
metric feature could appear in the spectra. The conduc- 
tance fluctuation is most enhanced for the contact show- 
ing a conductance value around 0.5 Gn, if the conduc- 
tance channel is a single one channel [2j. Therefore, the 
"non-symmetric" spectra were most frequently observed 
for the Au contact in H2 and D2 environment showing 
conductance of 0.1-0.5 Go- 

In order to determine the vibration energy from ac- 
tion spectra, 122 differential conductance spectra show- 
ing peaks in dl/dV spectra were collected for junctions 
having conductance of 0.05-0.4 Go- FigureO^b) shows the 
distribution of vibration energy. The broad features were 
observed around 33 and 66 meV. The two modes were 
also observed for the "normal" dl/dV spectra. The vi- 
bration mode of 33 meV was observed in Fig.[2](a). Since 
the phonon energy of the Au atomic contact is around 
10-20 meV 0], the observed vibration modes would be 
the vibration modes of the bonding between Au and hy- 
drogen. The distribution of the vibration energy could 
not be defined for the Au contact in D2 environment due 
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to the difficulty in obtaining the "peak" spectra. 

The structure of the Au atomic contact in H2 envi- 
ronment is discussed based on the present experimen- 
tal results and previously reported theoretical calcula- 
tion result [n| ■ After introduction of H 2 , the Au contact 
showed a conductance value below 1 Gq. The vibration 
modes between Au and hydrogen were observed around 
33 and 66 meV in the action spectra. The interaction 
between Au atomic contact with H atom or H2 molecule 
have been investigated by DFT calculation [l9( . H atom 
and H2 molecule were stably incorporated into the Au 
atomic contact. Because of the incorporation of hydro- 
gen into the contact, the conductance of the Au atomic 
contact decreased from 1 Go to 0.6-0.01 Go, depending 
on the atomic configuration of the contact. The decrease 
in conductance was possibly due to scattering or inter- 
ference of conducting electrons in the contact. No prefer- 
ential atomic configurations were found for the hydrogen 
incorporated Au atomic contact. Although there was no 
clear difference in conductance between the Au atomic 
contact with a single H atom and H2 molecule, there was 
clear difference in the vibration energy between two con- 
tacts • The Au atomic contact with a H atom showed 
the vibration mode around 150-220meV corresponding to 
the vibration of the H atom along the contact axis. Com- 
paratively, the Au atomic contact with a H2 molecule 
showed two transverse modes of the H2 molecule around 
25-150 meV, in addition to the internal H2 stretching 
mode around 180-250meV. These experimental and the- 
oretical calculation results indicated that a H2 molecule 
bridged Au electrodes in the present experimental con- 
dition. The formation of the Au/B^/Au junction was 
supported by the conductance trace. In the return con- 
ductance trace, the conductance jumped to a conduc- 
tance value of 0.1-0.4 Go- Before making the contact, H 2 
molecules would adsorb on Au electrodes without dis- 
sociation. In making the contact, a single B2 molecule 
would probably bridge Au electrodes in the initial stage. 
The vibration mode was also observed around 30-50meV 
for the contact having conductance of 0.1-0.4 Go, which 
was formed in making the contact. The close agreement 
in conductance and vibration energy supported that the 
atomic contact (conductance: 0.1-0.4 Go, vibration en- 
ergy: 40 meV) formed in breaking the contact would be 
similar to that in making the contact, and that a single 
H2 molecule would bridge Au electrodes. Now, the Au 
atomic contact in B2 environment was well character- 
ized by the conductance measurement, action spectra of 
the single molecule junction, and theoretical calculation. 
Briefly, the effect of surrounding hydrogen molecules on 
the Au/H2/Au junction is discussed. In the present ex- 
perimental condition, hydrogen molecules would adsorb 
on the Au/B2/Au junction. Since the interaction be- 
tween the hydrogen molecule and the hydrogen molecule 
in the Au/B2/Au junction would be small, the conduc- 
tance and vibration energy of the Au/B2/Au junction 



would be insensitive to the presence of the surrounding 
hydrogen. 

Next, dl/dV spectra of the Au/H 2 /Au and Au/D 2 /Au 
junctions were discussed. Three types of spectra ("nor- 
mal", "peak" and "non symmetric") were observed in 
dl/dV spectra. There was clear difference of the distri- 
bution of spectra between the Au/B 2 /Au and Au/D 2 /Au 
junctions. The percentage of the "normal", "peak" and 
"non symmetric" spectra were 7 %. 20 %, 73 % for 
the Au/B 2 /Au junction, and 6 %, < 1%, and 94 % 
for the Au/D 2 /Au junction, respectively. The "peak" 
spectra (action spectra) were frequently observed for the 
Au/B2/Au junction, while they were not observed for 
the Au/D 2 /Au junction. 

The peak in the dl/dV spectra is discussed using 
a model of vibrationlly induced two level systems (see 
Fig. [4]) [ll|]. In this model, the potential curve of the 
molecular junction is represented as the double well po- 
tential with ground states \&i and meta stable ^2 in the 
two energy minima. The two energy minima are sepa- 
rated by the activation barrier (Eac)- The molecular 
junction could be vibrationally excited by the conduc- 
tion electron. If the junction is fully exited, the junction 
with ground state ^1 could overcome the activation bar- 
rier, and change into the junction with meta stable ^2, 
which leads to the abrupt change in conductance (peak in 
the dl/dV spectra). The ^1 and ^2 states have slightly 
different local geometrical configurations, such as adsorp- 
tion site of a H2 molecule to Au electrodes, molecule tilt 
angle, configuration of the Au electrodes. Appearance 
of the "peak" spectra only for the Au/H^/Au junction 
indicated that the structural change could occur only for 
the Au/H 2 /Au junction. 

The difference between the two single molecule junc- 
tions is discussed by considering the zero point energy 
of the single molecular junction. Since a B2 molecule 
has half the mass of a D2 molecule, the energy of AU-B2 
vibration mode and zero point motion would be larger 
than those for D2. On the other hand, the potential 
curve would be the same for both single molecular junc- 
tions 15]. Therefore, Eac would be smaller for the 
Au/B 2 /Au junction than that for the Au/D 2 /Au junc- 
tion. The structural change could, thus, easily occur for 
the Au/H 2 /Au junction due to smaller Eac- The differ- 
ence in the dl/dV spectra between the Au/B2/Au and 
Au/D2/Au junction could be also explained by the num- 
ber of the excitation process accompany the structural 
change. If Eac was larger than the energy of the Au- 
H2 (D2) vibration mode, the system should be excited 
into the states with higher vibrational quanta in a ladder 
climbing manner in order to overcome activation barrier 
for the structural change [lij]. In the multiple excitation 
process, the excitation of the system would be getting 
harder with the number of the excitation process. The 
number of the excitation process to overcome the activa- 
tion barrier would be smaller for the Au/H2/Au junction 
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than that for the Au/D2/Au junction due to the larger 
energy of the AU-H2 vibration mode. Therefore, the 
structural change could easily occur for the Au/H 2 /Au 
junction. The above discussion was supported by the 
previously reported result Thijssen et al. reported 
that the "peak" spectra were observed for the Au con- 
tact in both H2 and D2 environment. They observed the 
"peak" spectra for the junction having a conductance 
value around 1 Go- In the present study, the dl/dV 
spectra were taken for the junction having a conductance 
value below 1 Go (around 0.1-0.4 Go). It is because the 
clear difference in the conductance histogram between 
clean Au contact and Au contact in hydrogen environ- 
ment was observed in the conductance regime around 
0.1-0.4 Go, and thus, the interaction between Au atomic 
contact and hydrogen was clear for the contact show- 
ing conductance around 0.1-0.4 Go- The "peak" spectra 
could not be observed for the Au/D2/Au junction. In the 
case of the vibrational heating mechanism, the reaction 
rate (R) per electron can be represented as R oc I n , where 
I and n are the current and the order of the reaction [3] . 
The reaction rate shows nonlinear power-law dependence 
on current. The structural change would occur for the 
Au/D 2 /Au junction having a high conductance value, 
while it would not occur for those having a low conduc- 
tance value. These experimental results suggested that 
the structural change in the Au/D2/Au junction would 
be the multiple excitation process. 



CONCLUSIONS 

The formation of the single Au/H^/Au and Au/D2/Au 
junction was revealed by the conductance measurement 
and action spectroscopy. In the dl/dV spectra, sym- 
metric peaks were observed for the Au/H2/Au junc- 
tion, while symmetric peaks were not observed for the 
Au/D2/Au junction. The peaks in the dl/dV spectra 
originated from the structural change induced by the 
phonon excitation. In the case of the Au/H2/Au junc- 
tion, the structural change induced by the phonon exci- 
tation could easily occur because the activation energy 
for the structural change and the number of excitation 
process would be small due to the large zero point energy 
and large energy of the AU-H2 vibration mode. 
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FIG. 1: (a) Conductance trace and (b) conductance histogram 
of Au contacts before (dotted line) and after (line) introduc- 
tion of H2 at the bias voltage of 100 mV. 



1 



0.15 

0.12 

0.09 

0.06 

0.03 
0.003 

0.000 

■0.003 




1 in H2 in D2 

7% 6% 



7i) 




-60 -30 30 

Bias Voltage (mV) 



60 




-60 -30 30 60 

Bias voltage (mV) 



FIG. 2: Schematic dl /dV spectra of Au contacts in H2 envi- 
ronment: (a) normal, (b) peak, (c) no symmetric feature type 
spectra. The percentage of the normal, peak, non symmet- 
ric type dl/dV spectra for single Au/H2/Au and Au/D2/Au 
junctions is shown in the figure. The number of the total 
dl/dV spectra was 200. 
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FIG. 3: Distribution of phonon energy for the Au/Hh/Au 
junctions obtained from the action spectroscopy of the single 
molecule junction. 
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FIG. 4: Potential energy surface of Au/H2/Au and 
Au/D2/Au junctions. 



